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Abstract. We consider a class of second order ordinary differential equations de- 
scribing one- dimensional systems with a quasi-periodic analytic forcing term and 
in the presence of damping. As a physical application one can think of a resistor- 
inductor-varactor circuit with a periodic (or quasi-periodic) forcing function, even if 
the range of applicability of the theory is much wider. In the limit of large damp- 
ing we look for quasi-periodic solutions which have the same frequency vector of the 
forcing term, and we study their analyticity properties in the inverse of the damping 
coefficient. We find that already the case of periodic forcing terms is non-trivial, 
as the solution is not analytic in a neighbourhood of the origin: it turns out to be 
Borel-summable. In the case of quasi-periodic forcing terms we need Renormalization 
Group techniques in order to control the small divisors arising in the perturbation se- 
ries. We show the existence of a summation criterion of the series in this case also, 
but, however, this can not be interpreted as Borel summability. 



1. Introduction 

Consider the ordinary differential equation 

ex + x + ex 2 — ef(ut), (1.1) 
where u) £ M. d is the frequency vector, f{tp) is an analytic function, 

/(</>) = ]T c^U, (1.2) 

with average a — a 2 , with a > (hence (/) = /o = a), and e > is a real parameter. Here and 
henceforth we denote with ■ the scalar product in M. d . By the analyticity assumption on / there are 
two strictly positive constants F and £ such that one has |/^| < Fe - ^"! for all u £ 1 . 
By writing 7 = 1/e the equation becomes 

x + jx + x 2 = f(u>t), (1.3) 

which describes a nonlinear electronic circuit, known as resistor-inductor-varactor circuit, subject 
to a quasi-periodic forcing function. Taking d = 1 and f(ujt) — a + /3sini, this equation has been 
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studied in Ref. [5], where, among other things, it has been found numerically that for 7 large enough 
there exists only one attracting periodic orbit and the corresponding period is 2it/uj = 2tt, the same 
of the forcing term. Furthermore one can prove analytically that such a periodic orbit is the only 
one in a neighbourhood of radius 0(1/7) around the point (a, 0). 

Here we give some further analytical support to such numerical findings. In particular we show 
that, if we take as forcing term an analytic periodic function, 

/(V) = ^e^/„ fo = a>0, (1.4) 

then for e small enough there is a 27r/cj-periodic solution, but this is not analytic in e — I/7 in a 
neighbourhood the origin in the complex e-plane. We find that such a solution is Borel-summable. 

We also show that by considering quasi-periodic forcing terms, as in (1.3), we still have a quasi- 
periodic solution with the same frequency vector u> of the forcing term, but we can only say in 
general that such a solution is analytic in a domain with boundary crossing the origin. 

Finally we shall see that considering more general nonlinearities introduces no further difficulties, 
and equations like 

x + jx + g{x) = f{ut), lim M^l = 00, (1.5) 

\x\^oo \x\ 

with g and / both analytic in their arguments, can be dealt with essentially in the same way. Simply, 
we have to impose a non-degeneracy condition on the function g, which reads as 

3x such that g(xo) = fo and g'(xo) 7^ 0. (1.6) 

In the particular case of homogeneous g(x), that is g(x) = ax p , with p > 2 an integer and ireK, 
the condition is automatically satisfied if p is odd (for any value of a) , while it requires afo > for 
p even, - as assumed in (1.1). 

The paper is organized as follows. For expository clearness we start with the case of periodic 
forcing terms. In Sections 2 and 3 we show that a periodic solution with frequency uj in the form of 
a formal power series in e (perturbation scries) is well defined to all orders, and it admits a natural 
graphical representation. In Section 4 we study further such a series, and we see that there is strong 
evidence to expect that it diverges. The best bounds that we are able to provide for the coefficients 
grow as factorials. To obtain bounds which allow summability of the perturbation series we have 
to perform a suitable summation in order to give the series a meaning. This is done in Section 5, 
and the resummed series is found to represent a 27r/o>-periodic solution which is Borel-summable in 
e. To prove the latter property we rely on Nevanlinna's improvement of Watson's theorem [12]. In 
Section 6 we consider the case of quasi-periodic forcing terms. We find that the perturbation series 
is well defined if the frequency vector of the forcing term satisfies a Diophantine condition, and, 
by using Rcnormalization Group techniques in order to deal with the small divisors problem, we 
find that the resummed scries still converges to a quasi-periodic solution, and it defines a function 
analytic in a domain containing the origin in its boundary. We shall see that the bounds we find 
do not allow us any more to obtain Borel summability, unlike the case of periodic forcing terms. In 
Section 7 we discuss how to extend the analysis to more general nonlinearities g(x), by requiring 
the condition (1.6) to be satisfied. 

The interest of the approach we propose is that it allows using perturbation theory which can be 
very natural in problems in which a small parameter appears. In fact analyticity in e for e close to 
(that is in 7 for 7 large enough) could be proved very likely with other techniques, but a naive 
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expansion in powers of e is prevented by the lack of analyticity in a neighbourhood of the origin. 
On the other hand the perturbation series gives a very accurate description of the solution, hence it 
is important to know that such a series is an asymptotic series, and its use is fully justified. Finally 
we can mention that the quasi-periodic solution we investigate is of physical relevance, hence it can 
be useful to study its properties. For instance in the case of the aforementioned resistor-inductor- 
varactor circuit in Rcf. [5] the 27r/w-periodic solution is numerically found to attract any trajectory 
which remains bounded in phase space. 



2. Formal analysis 

Consider first (1.1) for d = 1, that is 

ex + x + ex 2 =ef(ujt), (2.1) 

with /(f/0 given by (1.4). We look for bounded solutions (if any) which are analytic in e, that is of 
the form 

oo 

^) = E^ W W- ( 2 - 2 ) 

fc=0 

Inserting (2.2) into (2.1) and equating terms of the same Taylor order we find the set of recursive 
equations 

i(°> - 0, 

- _ T (0) _ T (0)2 , f 

jj — x dj i j i (2 3) 

iW = -jg(fc-i) _ J2 * {kl) x {k2 \ k>2. 

ki+k2—k—l 

From the first equation (zeroth order) we obtain that has to be constant, say x^ — c with Co 
to be determined. The second equation (first order) can give a bounded solution only if —c^ + a = 0, 
which fixes c = ^fa = a and gives x^ (t) as a periodic function with the same period of the forcing 
term: 

x« (t) = x« (0) + f dt' (/(wt 7 ) - a) . (2.4) 
Jo 

As each x^ (t) depends on the functions x^ k ) (t) with k' < k, we expect that if there is any periodic 
solution then it must have the same period as the forcing term. 

To continue the analysis to all orders it is more convenient to write the recursive equations (2.3) 
in Fourier space. The analysis to first order and the considerations above motivate us to write in 
(2.2) 

OO OO 

x(t) = £kx(k) (*) = £k e<I/U "4 fe) . (2-5) 

k=0 k=0 v£Z 

which inserted into (2.3) gives for v ^ 

4 0) = o, 

x<U - — 

iujv (2. 6 ) 

fcl+/c 2 — k — 1 ^1+^2— V 
k l, k 2 > D 
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provided that one has for v = 



= -4°> 2 + / , 

°= E E ^ 1} ^ 2 \ *>i. ( 2 - 7 ) 

fel+fc 2 =fe Vl+l>2=0 
k 1 ,fc 2 >0 



If we set = Ck then the first of (2.7) fixes, as already noted, 

c = a = y/a, (2.8) 
because one has /o = a > 0, while the second of (2.7) gives 

e e -r fe v fe :i = o. (2.9) 

fe'=o 

The latter equation, by taking into account (2.8) and the first of (2.6), can be more conveniently 
written as 



ci = o, c * = -fEE i 2"*' )it 2> ^ 2 < ( 2J °) 

c ° k'=i„ ie z 

which provides an iterative definition of the coefficients Ck as the right hand side depends only on 
the coefficients cy with k' < k. To deduce c\ — we used the first of (2.6), which, inserted into 
(2.9) for k = 1, gives 2coCi = 0, hence c\ = as c ^ 0. 
The following result holds. 

Lemma 1. Consider (2.1) with f given by (1-4) Then there exists a formal power series solution 
(2.2) whose coefficients x^ k \t) are analytic in t. If f is a trigonometric polynomial, that is in 
(1-4) one has \v\ < N for some N € M, then for all k > the functions x^{t) are trigonometric 
polynomials of order [(k + 1)/2}N, where [■] denotes the integer part. This means that one has 
xd k) = and x { ^ 1) = for \v\ > kN. 

Proof. The existence of a formal solution (2.2), with coefficients x^ k \t) analytic in t for all k > 0, 
follows from the analysis above. If / is a trigonometric polynomial of degree N, that the coefficients 
x^ arc trigonometric polynomials with the stated properties can be proved from (2.6) by induction 
on k. ■ 

Then the functions x^ (t) are well defined to all orders. Before discussing the issue of convergence 
of the formal power series defining such functions we look for a graphical representation of the 
coefficients 



3. Graphical representation and tree formalism 

We start by giving some abstract definitions. 

Definition 1 (Trees). A tree is a graph, that is a connected set of points and lines, with no 
cycle, such that all the lines are oriented toward a unique point which has only one incident line. 
Such a point is called the root of the tree. All the points in a tree except the root are denoted nodes. 
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The line entering the root is called the root line. The orientation of the lines in a tree induces a 
partial ordering relation between the nodes. We denote this relation by -<: given two nodes v and 
w, we shall write w <v every time v is along the path (of lines) which connects w to the root. 

Given a tree 9, we can identify the following subsets in 9. 

Definition 2 (Endpoints). We call E{9) the set o/endpoints in 0, that is the nodes which have no 
entering line. The endpoints can be represented either as white bullets or as black bullets. We call 
E\y{9) the set of white bullets and Eb{9) the set of black bullets. Of course Ew{9) U Eb{9) = E{9). 
With each v G Ew(9) we associate a mode label v v — ; an order label k v G Z + and a node factor 
F v = Ck v ■ With each v G Eb{9) we associate a mode label v v G Z\ {0}, and a node factor F v = f Vv . 

Definition 3 (Lines). We denote with L{6) the set of lines in 9. Each line £ G L{9) leaves a point 
v and enters another one which we shall denote by v' . Since £ is uniquely identified with v (the point 
which £ leaves), we may write £ — £ v . With each line £ we associate a momentum label i/<eZ and 
a propagator 

t lj vi = 0, 

and we say that the momentum vt flows through the line £. The modes and the momenta are related 
as follows: if £ = £ v one has 



i=l raeE B (fl) 
W-iv 



where £\, . . . , £ Sv are the lines entering v. 

Definition 4 (Vertices). We denote by V{9) the set of vertices in 9, that is the set of points which 
have at least one entering line. IfV{9) ^ we call the vertex v connected to the root the last vertex 
of the tree. If s v denotes the number of lines entering v call max^gy^) s v the branching number. 
One can have either s v = 1 or s v = 2. We set V s (9) = {v G V(9) : s v — s} for s = 1, 2; of course 
Vi(0) U V 2 {9) = V(9). We define also V (9) = {v e V{9) : v K = 0}; one has V (6) C V 2 (6). We 
require that either Vo(9) = or Vo(9) = {v n }, and that one can have v G V\ (6) only if V( v ^ 0. We 
associate with each vertex v G V(9) a node factor 

■1, s v = 2 and v ^ Vo(^) , 

F v = { -l/2c , s v = 2 and v G V (9) , (3.3) 
-{iwvz v ) 2 , Sy = I, 

which is always well defined as c ^ 0. 

We call equivalent two trees which can be transformed into each other by continuously deforming 
the lines in such a way that they do not cross each other. 
Let Tk tV be the set of inequivalent trees 9 such that 

(1) the number of vertices, the number of black bullets and the order labels of the white bullets are 
such that we have 

( fci + fc 2 + fc 3 = k, if v ^ , , , 

\k 1 + k 2 + k 3 = k + l, ifi/ = , 1 > 

if we set fa = \V(9)\, k 2 = \E B (9)\ and fc 3 = Y,veE w (d) k v, and 

(2) the momentum flowing through the root line is v. 
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We refer to Tk M as the set of trees of order k and total momentum v. 
With the above definitions the following result holds. 

Lemma 2. For all k > 1 and all v ^ one has 

4 fe) - £ Val(0), Val(0)=( J] 9t) f II F " ) ' ( 3 ' 5 ) 

oer ti „ \£eL(e) / \veE(e)uv(e) 

where Val : 7fc ; „ — > C is called the value o/ i/ie £ree. For fc > 2 and = one /las 

4 fe) = Cfc = *y^no), (3.6) 

eer fc , 

w/iere * means that there are two lines entering the last vertex vo of 9 and neither one exits from 
an end-point v with order label k v = 0. 

(k) (l) 
Proof. We can represent graphically Xq ' = Ck as in Figure 3.1a, xl, ', v ^ 0, as in Figure 3.1b, and, 

more generally, x^ as in Figure 3.1c. 

(fc) (1) (*) 

(a) (b) # (c)- 



Figure 3.1. Graphical representation of , x\, and xj> . For v = the latter reduces to 
the first graph, while for fc = 1 and v ^ it reduces to the second graph. In the first graph 
the momentum is not shown as it is necessarily v = 0. 



Then the third equation in (2.6) can be represented graphically as in Figure 3.2, if we associate 
with the nodes and to the lines the node factors and the propagators, respectively, according to the 
definitions (3.1) and (3.3). 




Figure 3.2. Graphical representation of the third equation in (2.6) expressing the coefficient 
for fc > 2 and v ^ in terms of the coefficients x^, ' with fc' < fc. In the last graph one 
has the constraints fci + fci2 = fc — 1 and v\ + V2 = v. 



Analogously (2.10) is represented graphically as in Figure 3.3, again if we use the graphical rep- 
resentations in Figure 3.1 and associate with the lines and vertices the propagators (3.1) and the 
node factors (3.3), respectively. 

(k) 

Note that in this way we represent graphically each coefficient xl in terms of other coefficients 
x v , , with fc < fc, so that we can apply itcratively the graphical representation in Figure 3.2 until 
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Figure 3.3. Graphical representation of the equation (2.10) expressing the coefficient c k for 
k > 2 in terms of the coefficients x^J ' with k' < k. Both k\ and k'2 arc strictly positive and 
ki + k2 = k\ moreover u\ + U2 = 0. 

only trees whose endpoints represent either XiP with v ^ (black bullets) or are left (white 
bullets). This corresponds exactly to the expressions in (3.5) and (3.6). ■ 

To get familiar with the graphic representation (3.5) and (3.6) one should try to draw the trees 
which correspond to the first orders, and check that the sum of the values obtained with the graphical 
rules listed above gives exactly the same analytical expression which can be deduced directly from 
(2.6) and (2.10). 

(2) 

For instance for k = 2 we obtain for x v , v ^ 0, the graphical representation in Figure 3.4 and 
for c 2 = .Tg 2 ^ the graphical representation in Figure 3.5. 



(0) 




Figure 3.4. Graphical representation of x v for v ^ 0. The second contribution has to be 
counted twice, because there is also a tree with the white and black bullets exchanged, - of 
course the latter has the same value. 




Figure 3.5. Graphical representation of C2 = Xq . There is no contribution with any white 
bullet carrying order label k = and k = 1 because of the restriction in the sum appearing 
in (3.6) and of the fact that ci = 0, respectively. 

For k — 3 we obtain for v ^ 0, the graphical representation in Figure 3.6 and for c 3 = x^ 
the graphical representation in Figure 3.7, where we have explicitly used that c\ = 0. 

And so on one can continue to higher orders. In general a tree 8 £ 7fe ;I/ looks like in Figure 3.8, 
where for simplicity no labels have been drawn other than the order labels of the white bullets. Note 
that each node can have only one or two entering lines, while the endpoints have no entering line at 
all. Moreover the momentum flowing through the line exiting a vertex v is equal to the sum of the 
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(0) 




(0) (0) (0) 




Figure 3.6. Graphical representation of x v for v ^ 0. The second and fourth contributions 
have to be counted twice, while the third one has to be counted four times. There is no 
contribution with any white bullet carrying the order label k = 1 as c\ = 0. 



(1) (0) (i) 




Figure 3.7. Graphical representation of C3 = Xq . The second contribution has to be counted 
twice, while the first one has to be counted four times. There is no contribution with any 
white bullet carrying the order label k = 1 as c\ = 0. 



momenta flowing through the lines entering v, according to (3.2): this is a sort of conservation law. 
The order of the tree is given by the number of vertices and black bullets plus the sum of the order 
labels of the white bullets minus the number of vertices in Vq{&). The latter is just |Vo(0)| = if 
G Tk,vi v + 0, and |Vu(0)| = 1 if 9 e T fc , . 

If a vertex v has s v = 1, that is it has only one entering line £, the latter can not come out from 
a white bullet. Indeed if this occurs one should have vg v = i/g = 0, hence F v = by (3.3), so that 
the value of the tree containing such a vertex is zero. 

Given a tree as in Figure 3.8 we can represent each white bullet according to the graphic rep- 
resentation in Figure 3.3, corresponding to the analytic formula (3.6), and expand again the two 
contributions and as sums of trees, and so on, iteratively, until the only white bullets 

which are left are the ones with order label k = 0. In this way we obtain a new graphic representa- 
tion where the trees still look like in Figure 3.8, but now there are a few differences: 

(1) all the white bullets v e Ew{0) have order labels k v — 0, and 

(2) there can be lines i € L{9) with momentum vi = which come out from vertices, that is Vq(0) 
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Figure 3.8. Example of tree appearing in the graphical expansions (3.5) and (3.6). The 
number of lines entering any vertex v can be only either s v = 1 or s v = 2, while no line enters 
the endpoints. The order of the tree is given by the number of elements in V{0) \ Vo(0) plus 
the number of elements in Eb(9) plus the sum of the order labels of the white bullets. Then, 
if fei, &2 and kz are the order labels of the white bullets in the figure the order of the tree is 
k = ki + k 2 + k 3 + 10 if Do <£ Vq(6) and k = ki + k 2 + k 3 + 9 if Do £ Vo(0). In the latter case 
one must have k\ > because of the constraint in the sum appearing in (3.6). 

can contain more than one or none element. 

Note that only lines coming out either from nodes in Vq{9) C V^(^) or from white bullets have 
vanishing momentum. 

The order of the tree is then given by the number of elements of V{9) U Eb{9) minus the number 
of elements of V {6), that is k = \V(9)\ + \E B {9)\ - \V (9)\. Of course v e V Q (9) if and only if the 
momentum of the root line is vanishing, that is 9 e 7^ for some k > 2. It is important to stress 
that no line entering a vertex v £ Vq (9) can come out from a white bullet (which now has necessarily 
an order label 0), because this would be against the constraint in the sum (3.6). This means that if 
two lines carrying zero momentum enter the same vertex v (so that v € Vq(9) according to (3.2)), 
then none of them can exit a white bullet. 

But up to these minor differences a tree representation like in (3.5) and (3.6) still holds. The 
advantage of these modified rules is that now the tree values are expressed no longer in terms of 
constants Ck to be determined, but only in terms of c which is known. A tree drawn according 
these new rules is represented as in Figure 3.8 with fci = k^ = k$ = (and in particular a tree of 
this kind can contribute only to with v ^ 0). Note that we could avoid drawing the order labels 
associated with the endpoints, as they are uniquely determined as k = for the white bullets and 
k = 1 for the black bullets. Of course, with respect to the caption of that Figure, now the order k is 
given by the number of elements in V(9) plus the number of elements in Eb{9) minus the number 
of elements in Vq{9). 

4. Formal solutions 

The sum over the trees in (3.5) and (3.6), with the new definition of the set 7^ .„ given at the end of 
Section 3, can be performed by summing over all possible "tree shapes" (that is trees without labels 
or unlabelled trees), and, for a fixed shape, over all possible assignments of mode labels. In the case 
of a trigonometric polynomial of degree N the latter can be bounded by (2N)\ E ( 6>> \, because each 
endpoint v can have either a mode label v v ^ 0, with \ v v \ < N, or the mode label v v — 0, while the 
case of analytic functions (or even to obtain bounds which are uniform in N) has to be discussed a 
little more carefully. The number of unlabelled trees with P nodes (vertices and endpoints) can be 
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bounded by 2 2P . 

Recall that V s {9) denotes the set of vertices v such that s v — s; of course V\{6) U V 2 {9) = V(6), 
and Vq{9) C V2{0). Analogously we can set 

L Q (6) = {£ E L(9) : n e = 0,} , 

Li(0) = {£ € L(6) :£ = £ v ,vG V^O)} , (4.1) 
L 2 (e)=L(e)\(Lo(9)UL 1 (9)), 

with the splitting made in such a way that one has 

n f "\\ n n w> \ n ^ n rL. 

«€Vi(e) teL^e) eeL 1 (8) eeL 2 {0) eeL 2 (8) 1 1 

/ i \ |vb(e)| 

n H^ur) ■ ! n ^|<^ w| . (4-2) 

veV o (0) ^ 0/ tiGB w (0) 

v£E B (9) vEE B (8) 

where for each line £ one has |^| < ^2 ve E w (0) \ u v\- 
The following result is useful when looking for bounds on the tree values. 

Lemma 3. Given a tree 9 with branching number s one has \E{9)\ < (s — 1)|V(#)| + 1. If k denotes 
the order of the tree9, that is \V(6)\-\V (6)\ + \E B (6)\ = k, one has the identity \L 1 (9)\ + \L 2 (0)\ = k, 
and the bounds | VI (6>) | < k, \V {6)\ < k-1, \E{9)\ < k and \E(9)\ + \V{6)\ < 2k - 1. 

Proof. It is a standard result on trees that one has J2vev(e)( s v ~ 1) = \E(0)\ — 1j so that the first 
bound follows. The bounds on |Vi(6% \V (6)\, \E{9)\ and \E(9)\ + \V{9)\ can be easily proved by 
induction, while the identity |Li(#)| + 1-^2 (#)| = k follows from the observation that all lines in 
L\(9) and L 2 {9) come out either from vertices or from black bullets, and they have non- vanishing 
momentum. ■ 

Hence the number of lines in L\{6) is bounded by k, so that in (4.2) we can bound 

( n n fc-^)<( n fb-w)( n ^ M/2k \^\) 

<€ii(e) veE B (8) veE B (9) leL^e) 

/9M IN* ( 43 ) 

<( J] Fe-^)f 2 -^ 

veE B (8) ^ ^ 

and in the second line the product can be used to perform the sum over the Fourier labels - and this 
gives a factor F k B\, with B 2 = 2e~£/ 2 (l - e - ^/ 2 ) -1 , - while the last factor is bounded by A\B k k\, 
for some constants A\ and B\. 

We can bound the value of a tree 9 by using the bounds (4.2) and (4.3), and Lemma 3. If we 
define 

Ei 1 = max{Bi, l^p 1 } max{c , FB 2 } max{l, (2c ) -1 }, (4.4) 

with Co = \foL, and take into account that the number of unlabclled trees in Tk, v is bounded by 
2 2fe ~ 1 (because each tree in 7fc ; „ has at most 2k — 1 nodes), then 
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,(*0 



< A^ k k\, 



,(*0 



(*) 



(4.5) 



where we have set £2 = £i2~ 2 . 

A bound like (4.5) is obtained also in the case of forcing terms which are trigonometric polynomials, 
because in general we can bound the factors \uve\ in (4.2) only with kN (see Lemma 1), and this 
produces an overall bound proportional to k\ Note that in that case the bound B2, arising from the 
sum over the Fourier labels, can be replaced with a factor 2N, and Bi can be replaced with |w|iV. 

Then we have proved the following result. 

Proposition 1. Given the equation (2.1) with f as in (1-4), there is only one periodic solution 
in the form of a formal power series, and the corresponding period is the same period 2tt/uj of the 
forcing term. The coefficients of such a formal power series satisfy the bounds (4-5). 

One could ask if the factorials arising in the bounds arc only a technical problem, or else they 
are a symptom that the series really diverges. To order k one can easily provide examples of trees 
which grow like factorials; see for instance the tree represented in Figure 4.1, where there are k — 1 
vertices with only one entering line. Then the corresponding value is 



Val(0) = { iu ;vf^-\f v = (iwv) k ~ 2 f v , 



(4.6) 



which behaves as k\ for large k. Furthermore it is unlikely that there arc cancellations with the 
values of other trees because the value of any other tree 9 E Tk,v can be proportional at most to 
(icju) p , with p < k — 2 (strictly). Hence we expect that the coefficients ui k \ even if well defined to 
all orders, grow like factorials, so preventing the convergence of the series. 

(1) 



Figure 4.1. Example of tree whose value grows as a factorial. If k is the order of the tree 
(hence there are k — 1 vertices and 1 black bullet), then the value of the tree is given in (4.7). 



The lack of analyticity is further supported by the following fact. If we consider (2.1) without the 
quadratic term and with a = 0, that is 

ex + x = f{urt), / = 0, (4.7) 

in Fourier space, we find x = and iuov{l + ieuov)x u = f v for v ^ 0. Hence the equation is trivially 
solvable, and it gives 

= } ; ^ (4.8) 

Of course the solution x(t) of the linear equation is not analytic in e (in a neighbourhood of the 
origin) when / is an analytic function containing all the harmonics, as each point e — i /uv represents 
a singularity point for x(t), and such points accumulate to the origin as v — > 00. Then it is likely that 
also when the quadratic terms are taken into account the solution can not be analytic. Therefore 
giving a meaning to the perturbation series requires some more work, that we are going to discuss 
next. 
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An important remark is that for any k > 1 there is no tree whose value can be bounded worse 
than proportionally to a factorial, as the estimates (4.5) show: indeed they have been obtained by 
bounding separately the value of each single tree. This observation will play an important role in 
the forthcoming analysis. 



5. Periodic forcing terms 

To deal completely with the case of analytic functions and prove existence of the periodic solution, 
we have to modify the graphical expansion envisaged in the previous sections. 

Let us come back to the equation (2.1), and write it in Fourier space. For v ^ and denoting 
with x v the i/-th Fourier coefficient, we obtain 

e(itjjvfx v + iujvx v + e x Vl x V2 = ef v , (5.1) 
provided that for v = we have 

x Vl x„ 2 = 0. (5.2) 

^1+^2=0 

Let us rewrite (5.1) as 

e{iujv) 2 x v + ibjvx v + \ie ^ x Vl x V2 = fxef v , (5.3) 
and look for a solution x{t) which is analytic in /z, which suggests us to write 

oo 

x{t) = Y.v kx[k]( t)- ( 5 - 4 ) 

k=0 

Of course we want that at the end the value fi = 1 is inside the analyticity domain. Note also that 
now x^ , the coefficient to order k, has a different meaning with respect to the previous expansion 
(1.4) in powers of e, and for this reason with use a different symbol to denote it. We shall call 
resummed series the series (5.4), because the coefficients x^(t) depend on e, and are given by the 
sum of infinitely many terms of the formal series (2.5). 

Again for k = we have to take xL°' = for v ^ and fix cq = x^ = y/a, with a = /o- 

To order k > 1 (in jj) we obtain for v ^ 

iuv(l + ieuv)xW =efj k ,i-s ^ ]T x^x^, (5.5) 

fcl+fc2 = fe-l 1>1+V2=V 

while for v = we require 

E E « 2] =0. (5.6) 

ki+k 2 =k vi+v 2 =v 

By setting = x\^ the latter equation can be written as (cf. (2.10)) 

ci = o, ^xEE 1 ^ 1 - 1 - ( 5J ) 
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Then we can proceed as in Section 2, with some slight changes that we now explain. First of all 
note that (5.5) gives for v ^ 



0, 

iojv(l + ieuv) ' (5.8) 

TTT~- s E E W/ 1 , 

iluv 1 + leujv) ^— ' ^ 1 2 

fcl+fc 2 = fc— 1 I>l+I>2=»' 

Then the graphical representations of xff' , rrl 1 ' and x$ arc as in the previous case, with the only 
change in the representation of the order labels (because of the square brackets instead of the 
parentheses); see Figure 5.1. 




Figure 5.1. Graphical representation of x l , x\,' and x\,\ For v = the latter reduces to 
the first graph, while for k = 1 and v ^ it reduces to the second graph. In the first graph 
the momentum is not showed as it is necessarily v = 0. 



40] 



On the contrary the graphical representation of the third equation in (5.8) is as in Figure 5.2. 




Figure 5.2. Graphical representation of the second equation in (5.8) expressing the coefficient 

\k\ \k'\ 
x v 1 for k > 2 and v ^ in terms of the coefficients x , with k < k. In the last graph one 

has the constraints fci + k2 = k — 1 and v\ + V2 = v. 



At the end we obtain a tree expansion where the trees differ from the previous ones as they contain 
no vertex with only one entering line. With the previous notations this means that L\(6) = and 
Vi(9) = 0, hence V{9) = V^(^). Moreover also the propagators and the node factors of the vertices 
are different, as (3.1) and (3.2) have to be replaced with 



and, respectively, 



1 1, vi = 0, 



v ~ 1 -l/2c , « G V (9) , l5 - iUj 



and we recall once more that only vertices v with s v = 2 are allowed. Finally, the node factors 
associated to the endpoints are F v — Ck v if v is a white bullet and F v — ef Uv if v is a black bullet. 
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Figure 5.3. Example of tree appearing in the new graphical expansion. The number of lines 
entering any vertex v can be only s v = 2. The order of the tree is given by \B(8)\ — |Vb(0)|. 
All the white bullets have order labels 0, as well as all the black bullets carry a label (1); 
hence we can avoid drawing explicitly such labels. 



As in Section 3 we can envisage an expansion in which all white bullets v have k v = (simply 
by expanding iterativcly in trees the white bullets of higher order). A tree appearing in this new 
expansion is represented in Figure 5.3. 

With the notations (4.1), we obtain the bounds 



n 

vev(e)\v a (0) 



< 



<m |vw| , i n 

\Vo (0)1 

. I n * 

vEE B (e) vEE B (0) 



n 



1 



\u3vA 1 + ieojuA 



n *•* 

vev (e) 



\E W (9)\ 



(5.11) 



where we have used again the bound < Fe"^"l, for suitable (strictly) positive constants F and 
£, which follows from the analyticity assumption on /. 
For real e we can bound each propagator by 



1 1 

\9e\ < i 1 < r 

so that the value of any tree 9 € 7fc ; „ can be bounded by 



|Val(0)| < |e| fc |wr' £ (max{co,F})' £ (max{l,l/2co})' £ J] e~^, 

v£E B (8) 



(5.12) 



(5.13) 



where we have used again Lemma 3. If we write 



veEb(o) \veE B (e) 



(5.14) 



we can proceed as in ection 4: we use the last product to perform the sum over the Fourier labels, 
which gives a factor B%, whereas the sum over the unlabelled trees gives a factor 2 2fe ~ 1 . At the end 
we obtain 

x^(t) <pq\ (5.15) 
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where we have set [i^ 1 = 4 max{l, l/2co} max{FB2, co}|e|. Hence the radius of convergence 
yu 01 the series expansion (5.4) is bounded as ^ > f-2 = 0(l/|e|), so that for e small enough, say 
\e\ < £3 = (4 max{l, l/2co} max{FB2, Co}) -1 , the value /x = 1 is inside the analyticity domain. 
We can summarize the results found so far as follows. 

Theorem 1. Given the equation (2.1) with f analytic, there exists e > such that for all real 
s with \e\ < £0 there is only one periodic solution which admits a formal expansion in powers of 
e, and the corresponding period is the same period 2ir/u of the forcing term. An explicit bound is 
£0 > £3 = 0(u). 

Note that if uo is very large then very large values of s are allowed. 

We can investigate further the regularity properties in e of the periodic solution found in Theorem 
1, and see what happens for complex values of e. 
We need the following preliminary result (see Figure 5.4a for the region Cr). 



Ime 



(a)- 



Ree 



(b) 



Ime 














Ree 







Figure 5.4. Region Cr in the complex £-plane (a) and strip-like region of analyticity Sb of 
the Borcl transform (b). The region C_r is the union of two discs of radius R/2 and centers 
(±R/2,0). 



Lemma 4. Given ui > and < R < 1/Aw let Cr be the pair of discs Cr = {e : |Rce x | > R x }. 
For all e e Cr and all v £ Z \ {0} one has \iuv(l + ieuov)\ > ui/2. 

Proof. Write e = a+ib and x = lov, so that one has \iuv{l + ieuji>)\ = \x\y/(l — bx) 2 + (ax) 2 = F(x). 
If e eCr one has \a\ > b 2 /2R. Fix < A < 1. If |1 - bx\ < A then ^(1 - bx) 2 + (ax) 2 > \ax\ > 
b 2 \x\/2R> \b\(l-A)/2R, so that F(x) > (1-A) 2 /2R. If |1 — 6sc| > A then ^(1 - 6a;) 2 + (ax) 2 > A, 
hence F(x) > A\x\ > uoA. Then choose A = 1 - \fusR > 1/2; this gives F(x) > uj/2. ■ 

Now fix < R < R = £3 so small that \oj\R < 1/4, and consider the corresponding domain Cr. 
We can apply Lemma 4 and deduce that any propagator gi is bounded by \ge\ < 2/\u>\ for all e € Cr. 
This allows us to obtain the following result. 

Proposition 2. There exists R > small enough such that in the domain Cr one has the asymptotic 
expansion 

N-l 

x (t) = £ k x {k) (t) + 9i N (e), |SHjv(e)| < AB N N!\e\ N , (5.16) 

k=0 
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where the constants A and B are uniform in N and in e. 

Proof. Write x(t) as x(t) = xjv(t) where x^(t) is given by the sum of the first N — 1 orders 

of the formal power series expansion of the solution x(t) as in (5.16). For e € Cr the function (5.4) 
with \i = 1 is C°° in e, hence we can estimate 9^7v(e) with a bound on the JV-th derivative of x(t) 
in Cr, and this gives the bound in (5.16). ■ 

Of course the constants A and B in (5.16) are explicitly computable; in particular one finds 
B = 0(e^). 

Then we are under the assumptions where Nevanlinna's theorem [11] (see also Ref. [12]) can be 
applied, and we obtain that the function 

oc 

< 8(t;e) = J2-e k xW(t) (5.17) 

k=0 

converges for |e| < B (with B given in (5.16)) and has an analytic continuation to Sb — {e : 
dist(e,lR+) < B} (see Figure 5.4b), satisfying for some constant K the bound |Q3(i;e)| < Ke^ E ^ R 
uniformly in every Sb' with B 1 < B. The function x(t) can be represented as the absolutely 
convergent integral 

x(t) = - / e- s / ef B(t;s)ds (5.18) 
£ Jo 

for all e e Cr, and this property can be stated by saying that x(t) is Borel-summable (in e) and 
03 (t; e) is its Borcl transform [9]. This implies that the function given by the summation procedure 
described in Theorem 1 is unique. Therefore we have obtained the following result, which strengthens 
Theorem 1. 

Theorem 2. The solution given by Theorem 1 is Borel-summable at the origin. 

Note that Watson's theorem can not be invoked to obtain this result because the singularities are 
along the imaginary axis. 

In particular if f(u)t) = a + /3sini then there is a periodic solution x(t) = a + sflcost + 0(e 2 ), 
with a = yfot, which has period 2ir and moves around the fixed point (x, x) = (a, 0), and close to it 
within 0(e). No other periodic solution analytic in e can exist. 

We conclude this section with two remarks. The summation criterion envisaged in this Section is 
reminiscent of that used (in a more difficult situation) in Ref. [6] for hyperbolic lower-dimensional 
tori. However in that case we have not been able to prove Borel summability because to order k 
the bounds are like (k\) a for some a > 1. Neither extensions to Watson's theorem [9] analogous to 
Nevanlinna-Sokal's result (as those developed in Ref. [4]) can be used because the exponent a is too 
large. We shall find a very similar situation in next section. 

The lack of analyticity in e in a neighbourhood of the origin is due to the accumulation of singu- 
larity points along the imaginary axis in the complex e-plane (where the quantity 1 + ieuv vanishes 
for v e Z). The analyticity domain is tangential to the imaginary axis, and this allows us to apply 
Nevanlinna's theorem. We find that this situation has some analogies with a different problem, the 
analyticity properties of rescaled versions of some dynamical systems, such as Siegel's problem [3] 
- and its linearization as considered in Ref. [10] -, the standard map [1] and generalized standard 
maps [2], for complex rotation numbers tending to rational values in the complex plane. In those 
cases, however, only non-tangential limits could be considered. Of course the situation is slightly 



16 



more complicated there, because the set of accumulating singularity points is dense - and not only 
numerable as in the present case. 



6. Quasi-periodic forcing terms 

In the case of analytic quasi-periodic forcing terms, we shall assume a Diophantinc condition on the 
rotation vector w, that is 

> C \u\- T ViveZ d \{0}, (6.1) 

where \v\ = \u\i = \vi \ + . . . + \vd\, and Co and r are positive constants. We need r > d— 1 in order 
to have a non-void set of vectors satisfying the condition (6.1), and r > d — 1 in order to have a full 
measure set of such vectors. For simplicity (and without loss of generality) we can assume Co < 7/2, 
with 7 = min{l, |c|}, where c is a suitable constant to be fixed asc = — 2co, with Co = \/a. 
The equation of motion can be written in Fourier space as 

iu> ■ v (1 + ieuj ■ u) x v + e ^ x Vl x V2 = ef u , (6.2) 

V 1 +U 2 =l' 

and the formal expansion for a quasi-periodic solution with frequency vector u) reads as 



x{t) = J2 £kx(k) (*) = E £k E . (6-3) 



'- c 

fc=0 fe=0 



and to see that the coefficients x\, ' are well defined to all orders k > one can proceed as in 
Section 2, with no extra difficulty. In particular the Diophantinc condition (6.1) is sufficient to 
assure analyticity in t of the coefficients x^ k \t). 

Also the graphical representation can be worked out as in Section 3. The only difference is that 
now the propagators of the lines with non- vanishing momentum v>i, which is defined according to 
(3.2), with the vectors replacing the scalars, are given by ■ ug), the node factors associated 

with the vertices v with s v = 1 are given by F v = —(ico ■ ve v ) 2 , and the node factors associated 
to the black bullets v are given by F v = f Vv , with u v e if \ {0}. All the other notations remain 
unchanged. 

This yields that the propagators and the node factors can be bounded as in (4.2) and (4.3), with 
just a few differences of notation. More precisely one has 

n p v\\ n H- n hn. i n H< n k^-^i^r^ 

/ 1 \ |vb(e)| 

n ^Uc-) . ! n ^\<^ m , (6.4) 

H F V \<F\ E °W [] e-«l"»l, 

veE B (8) veE B (9) 

where the only bound which introduces a real difficulty with respect to the case of periodic forcing 
terms is the second one in the first line. Indeed it is the source of a small divisors problem, which 
can not be set only through the Diophantine condition (6.1). 
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To each order k we obtain for x^(t) a bound like AB k k\ maK ^ 1 ' T \ where the factor 1 arises 
from the propagators of the lines in Li{9) and the factor t from those of the lines in L 2 (6) in 
(6.3). The last assertion is easily proved by reasoning as in (4.3), with max{|u>| |i/^|,CjJ~ |^| r } < 
maxjC^ 1 , |a;|}|i/£| max ^ 1,T ^ replacing vg. In particular only for d = 2 and r = 1 we obtain the same 
bound proportional to k\ as in the case of periodic solution (of course with different constants A 
and B). Note that the vectors satisfying the Diophantinc condition (6.1) with r = 1 for d = 2 is 
of zero measure but everywhere dense. An example of vector of this kind is u> = (1,70), where 
7o = (V5 — l)/2 is the golden section. 

However, to deal with the problem of accumulation of small divisors and discuss the issue of 
convergence of the series, we need Renormalization Group techniques. The first step is just to 
introduce a multiscale decomposition of the propagators, and this leads naturally to the introduction 
of clusters and self-energy graphs into the trees. The discussion can be performed either as in Ref. [6] 
or as in Refs. [8] (and in Ref. [7]). We choose to follow Ref. [8], which is more similar to the present 
problem because the propagators are scalar quantities and not matrices. In any case, with respect to 
the quoted reference, we shall use a multiscale decomposition involving only the quantities \u> ■ V(\, 
that is without introducing any dependence on e in the compact support functions. Indeed this is 
more suitable to investigate the analyticity properties in e, and, as we shall see, we shall not need 
to exclude any real value of e in order to give a meaning to the resummed series, a situation more 
reminiscent of Ref. [6] than of Refs. [8]. 

In the following we confine ourselves to outline the main differences with respect to Ref. [8] . Let 
us introduce the functions ip n and \n, f° r n > 0, as in Ref. [8], Section 5. In particular ?/>„(|x|) 7^ 
implies \x\ > 2~( n+1 )C and Xn(M) 7^ implies \x\ < 2~ u Cq. Wc shall define recursively the 
renormalized propagators gf^ — g^ n \u> ■ vf, e) and the counterterms M\ n \u> ■ u; e) on scales n as 

1, M^(x;e) 

MM) 



i°\x;e 



gW(x;e 
M [n] {x-e 
MW(x;e 



0, 

M^(x-e) 



ix(l + iex) ' 

Xo(\x\) ■ ■ .Xn-l(\x\)tp„(\x\) 



E vt(x ]£ ), 



-- 1 Tesl a 



ix{l 
M [n ' 

OO 



+ iex)+M^-^(x;e) ' 

e) + xo(M) • • • Xn-i(\x\) X n(\x\)M^(x;s), 



(6.5) 



fe=i T es£ n 



where the set of renormalized self-energy graphs S^ n and the self-energy graphs Vt(x; s) are defined 
as in Ref. [8], Section 6. We have explicitly used that the first contribution to the self-energy graphs 
is of order k = 1 (sec Figure 6.1). Note that one has Xo(M) • ■ • Xn-i(M)Xn(M) = Xn(M) ; so that 
if g [n] (x:e) ^ then one has 2-( n+1 )C < \x\ < 
Then one defines for k > 1 



E Val (^ 

eer fc ,„ 



Ck 



E * Val (^ 

eeT k:0 



where the tree value is defined as 



Val(6>) 




n *•* 

v£E(e)UV(6) 



(6.6) 



(6.7) 
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(0) 




+ 



Figure 6.1. First orders contributions to the counterterm arising from self-energy graphs 
of order k = 1 and k = 3. The (dashed) external lines do not enter into the definition of 
self-energy graph, and they have been drawn only with the aim of helping visualizing the 
structure of the self-energy graph. 

and, as before, * means that there are two lines entering the last vertex w of 9q and neither one 
exits from an endpoint v with order label k v — 0. For k = 1 the second of (6.7) has to be interpreted 
as Ci = 0. 
Furthermore one has 



M^{x;e) = M^(0;e) + O(e 2 x), 

Af[°](0;e) = -2ec + M[ ol (0;e), M' ol (0;e) = 0(e 2 ) 



(6.8) 



and an easy computation shows (cf. Figure 6.1) that 



«f (M - «■! - ,|) ( „ . _ + 0{ e% (6.9) 

so that in fact one has M2 '(0;e) = 0(e 3 ). 
Moreover to higher scales one has MW(i;e) = mW(0;e)|O(e 3 i), with 

MN( 0;£) = -.3l^ ^-^n^-u\) - ■ +Q(e 4 ), (6.10) 

SO that each MW(0;e) is a higher order correction to (0; e) and it decays exponentially in n 
(because of the compact support functions). 
The following result holds. 

Lemma 5. Assume that the renormalized propagators up to scale n — 1 can be bounded as 



Am] 



< C~ 1 2 l3ni (6.11) 



for some positive constants C\ and (5. Then for all n' < n — 1 the number N n i{9) of lines on scale 
n' in 9 is bounded by 

N n ,(0)<K2- n '/ T \"v\, (6-12) 

vEE B (8) 

for some positive constant K. If\e\<eo, with e small enough, then for all n' < n one has 

|Ml n '](x;e)| < Dxlepe-^ 2 " 7 ', \d x M^'\x; s)\ < D^ef ''e'^ 2 " ' '\ (6.13) 
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for some C\ -independent positive constants D\ and L> 2 . Only the constant D\ depends on (3. The 
constant e can be written as e = CiC^^C^, with C 2 and C 3 two positive constant independent of 
(3 and C\ . 

Proof. The proof can be easily adapted from the proofs of Lemma 1 and Lemma 2 of Ref. [8] . ■ 

So we are left with the problem of proving that the renormalized propagators satisfy the bounds 
(6.11). To this end let us introduce the notation 



F(x) = F (x) + ci(e) e + c 2 (e, x)e 2 x, F (x) = ix(l + iex), 



(6.14) 



with x = u) ■ i> and the functions ci(e) and c 2 (e, x) such that c\(e) = c + c 3 (e) e, with c ^ 0, and 
the functions |c 2 (£,x)| and |c3 (e)J bounded by a constant c' uniformly (in e and x). Recall that 
7 = min{l, \c\} and Co < 7/2. 

Fix A e [0, 1]. Set B R (0) = {e e C : |e| < R} and V R , X = {e = a + ib E B XR (0) : \a\ > X\b\} (see 
Figure 6.2). The following result refines Lemma 4. 




Ree 




Ree 



Figure 6.2. Region T>b.,\ in the complex e-plane for A = tan-7r/6 (a) and for A = 1 (b). One 
can write A = tan <p, where ip is the angle between the imaginary axis and the line a = \b. 



Lemma 6. Given < R < l/4Co let Cr be defined as in Lemma 4- For all e € Cr and all x one 
has 1^0(^)1 > niin{Co, |a;|}/2, while for all e e T>R t \ one has \F (x)\ > A|a;|/2. 

Proof. Write e = a + ib, so that |F (a;)| = \x\y/(l - bxf + (ax) 2 . For e E C R set A = 1 - VOTR. If 
\x\ > Co, for 1 1 — bx\ < A one has |F (x)| > |aa; 2 | > b 2 x 2 /2R > C /2, while for |1 - bx\ > A one has 
\F (x)\ > A\x\ > \x\/2 > Co/2. If\x\ < C , for |l-6a;| < A one has |F (x)| > |aa; 2 | > C /2 > \x\/2, 
while for |1 - bx\ > A one has \F (x)\ > A\x\ > \x\/2. For e e V R . X set A = 1/2: one finds 
\F {x)\ > X\x\/2. U 

Then the following result holds. 

Lemma 7. Set x = uj v and assume \x\ < Co- Then if R is small enough one has |.F(a;)| > A7|x|/8 
for all e € T>r^\. 

Proof. Set Fi(x) = F (x) + ce and e = a + ib. Then Fi(x) = i(x + b(c - x 2 )) + a(c - x 2 ), and 
\F( x)\ > |F i(x)| - c'|e| 2 (l + |x|). If \x + b(c - x 2 )\ > \x\/2 and \bc\ > 4|x| one has |Fi(ar)| > 
\c\\/b 2 + a 2 /2 = \cs\/2, so that \F(x )\ > \ce \/4 > \cb\/A > \x\. If \x + b(c - x 2 )\ > \x\/2 and 
\bc\ < 4\x\ one has |Fi(x)| > ^ m&x{\/ x 2 + a 2 , \e\/4}/2, so that |F(a;)| > 7V 1 x 2 + a 2 /A > 7|x|/4. If 
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\x + b(c-x 2 )\ < \x\/2 one has |6(c-x 2 )| > \x\/2 and \bc\ < 3\x\, which give |e| 2 < 3\e\Va 2 +x 2 /~f < 
3XR(\a\ + \x\)/j, and |Fx(x)| > \a(c - x 2 )\ > \a(c - x 2 )\/2 + (X\x\/2)/2 > >y\(\a\ + |a;|)/4, so that 
\F(x)\ > 7 X\x\/8. m 



Then we can come back to the bounds of the renormalized propagators, and prove the following 
result. 

Lemma 8. If R is small enough for all n > and all e € I?r,a the renormalized propagators 
g^(x;e) satisfy the bounds (6.11) with (3=1 and C\ = XC4, with a X-independent constant C4. 

Proof. The proof can be done by induction on n. For n — the bound is trivially satisfied by 
Lemma 6. Assuming that the bounds hold for all n' < n then we can apply Lemma 5 and deduce 
the bounds (6.13). In turns this implies that the renormalized propagators on scale n can be written 
as g[™](x;e) = 1/F(x), with F(x) written as in (6.14) for c = — 2co (cf. (6.8)), and for suitable 
functions c\{e) and C2(s, x), depending on n and satisfying the properties listed after (6.14) for some 
n-independcnt constant c'. Then by Lemma 7 the renormalized propagators g^ n \x;s) satisfy the 
same bounds (6.11) with C\ = 0(A) for e e T>^\. ■ 

Of course for real e the bound (6.11) is trivially satisfied, with 0i = 2~ 1 O - This follows from 
Lemma 8 with A = 1, but it is obvious independent of that result because one has ci(e) = ce + 0(e 2 ), 
with c = — 2co € M.. If we want to take also complex values of e, we have analyticty in a domain D 
which can be written as T>r = Uxe[o,i]f.R,A- One can easily realize that the region Cr is contained 
inside the domain V R (cf. Figure 6.3). Fix A = tantp, with ip e [0, 7r/4] (see Figure 6.2): for all such 
ip the line which forms an angle ip with the imaginary axis (see Figure 6.2) and passes through the 
origin intersects the boundary of T>r at a distance R tan ip from the origin and the boundary of Cr at 
a distance R sin ip. Hence we have an analyticity domain of the same form as in the case of periodic 
forcing terms. Nevertheless the results found so far do not allow us to obtain Borel summability, 
notwithstanding a circular analyticity domain Cr is found, as the bounds which are satisfied inside 
the region Cr are not uniform in e (because of the dependence on A) . 



Note hat f3 = 1 in (6.11) is the same exponent appearing in the bounds of the propagators in the 
formal expansion. To obtain uniform bounds in a domain Cr, for some value of R, we have to allow 
larger values of (3. The following result is obtained. 

Lemma 9. Set x = uj ■ u and assume \x\ < Co- If R is small enough one has \F{x)\ > ^y\x\ 2 /2 for 



Ime 
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Figure 6.3. Regions T>r and Cr in the complex e-plane. The first one is the grey region, 
while the second one is the region contained inside the two circles. 



21 



all e S Cr. 



Proof. Set Fi(x) = F (x) + ce and e = a + ib. If \x + b(c - x 2 )\ < \x\/2 one has \b(c - x 2 )\ > 
\x\/2 and 3\x\ > \bc\ > \x\/4. Hence \a(c - x 2 )\ > b 2 \c - x 2 \/2R > |a;| 2 /16i?|c|, so that one has 
\Fi(x)\ > \a(c -x 2 )\ > \ac\/4 + \a(c - x 2 )\/2 > j(\a\ + x 2 /16Rc)/2. On the other hand one has 
\e\ 2 = a 2 + b 2 < a 2 + 9x 2 /c 2 , so that \F{x)\ > |iq(a;)| - 2c'\e\ 2 > |Fi(a;)|/2 > jx 2 /2. The case 
\x + b(c — x 2 )\ > \x\/2 can be discussed as in Lemma 7, and it gives > 7|x|/4. ■ 

Then we can prove the following result by proceeding exactly as in the proof of Lemma 8. 

Lemma 10. If R is small enough for alln> and all e £ Cr the renormalized propagators (x; e) 
satisfy the bounds (6.11) with {3 = 2 and C\ a suitable constant. 

The advantage of Lemma 8 with respect to Lemma 10 is that the bound of R is better, which 
means that the domain Cr contained inside Vr in the first case is larger than the domain Cr of 
the second case. The advantage of Lemma 10 is that it allows to obtain uniform bounds inside the 
corresponding domain Cr. Nevertheless, because of the factor [3 = 2, a bound AB k k\ 2r is obtained 
for the coefficients x^ (t) of the formal solution, and a result analogous to Proposition 2 can be 
proved also for the present case, with 7V! 2r replacing iV!; we do not give the details as the proof 
is identical. Hence the bounds that we have are not good enough to obtain Borel-summability in 
the case of quasi-periodic forcing terms, a situation strongly reminiscent of that encountered in 
Ref. [6]. In fact at best one can set t = 1 for d = 2 (which, as noted above, corresponds to a set 
of Diophantine vectors of zero measure but everywhere dense), but this in turn implies a bound 
proportional to AH 2 , which is not enough to apply Nevanlinna's theorem. 

The conclusion is that the resummed series 

oo 

x(t)=J2» k x lk] (t), (6-15) 

k=0 

where the coefficients x^ (t) are given by 

xW(t) e^*xl fcl , (6.16) 
ve1 d 

with x$ defined by (6.6), is well defined and converges. In general it is not obvious - even if 
expected, - that (6.15) solves the equation of motion (1.1). Indeed, unlike the case of periodic 
forcing terms, we have no result, such as Nevanlinna's theorem on Borel summability, which we can 
rely upon in order to link the resummed series to the formal series. Therefore we have to check 
by hand that by expanding in powers of e the resummed series we recover the formal power series 
(6.3). This means that the resummed scries, which in principle could be unrelated to the equation 
of motion (because of the way it has been defined), in fact solves such an equation. Such a property 
can be proved by reasoning as in Ref. [8], Section 8. Again we omit the details, which can be easily 
worked out. 

We can summarize our results in the following statement. 

Theorem 3. Given the equation (1.1) with f analytic in its argument and oj satisfying the Dio- 
phantine condition (6.1), there exists e such that for all real e with \e\ < £q there is a quasi-periodic 
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solution with the same frequency vector of the forcing term. Such a solution extends to a function 
analytic in a domain V R like in Figure 6.3, with R — e . 

The conclusion is that the summation criterion described here gives a well defined function, which 
is quasi-periodic and solves the equation of motion (f .1), but the criterion is not equivalent to Borel 
summability any more. In particular the issue if such quasi-periodic solutions are unique or not 
remains open, as in Ref. [6]. 



7. Extension to more general nonlinearities 

When considering the equation (1.5) the formal analysis of Section 2 (and of Section 6 in the case 
of quasi-periodic forcing terms) can be performed essentially in the same way. If we write 

[s(-)]t fc) = E^ E k>o, (7 ' 1} 

p=0 V' ki + ... + k p =k 
v-l+...+v p =ii 

then the recursive equations for v ^ arc 



= 0, 

(1) _ fu 

ILJ ■ V 

xl k) = -(t« • u)xi k - 1] - _!— fo^- 1 ), k > 2, 



(7.2) 



while the compatibility condition becomes [ff(a:)]o^ = /o<5fc,o for k > 0. The latter for k = 
gives g(co) = fo, while for k > 1 gives g'(co)ck + R{co, ci, . . . , Ck-i) = 0, where the function 
R(cq, ci, . . . , Cfc_i) depend on the coefficients to all orders k' < k, hence, in particular, on the 
constants c , . . . , Ck-i- Therefore the constants can be fixed iteratively as 

Ck = 77 — T- R ( c o> c i> • • • > c fc-i); (7-3) 

provided that one has g'(co) ^ 0, so that under the conditions (1.6) one has the formal solubility 
of the equations of motion (1.1). Note that the fisrt condition in (1.6) requires fo G Ran(g), and if 
such a condition is satisfied then the condition on the derivative is a genericity condition. Note also 
that the class of functions g(x) which are not allowed depends on / (more precisely on its average 
.fo)- For instance an explicit example of function which does not satisfy (1.6) is g(x) = 3x 2 — 2a; 3 if 
fo = I- 

The graphical representation differs from that of the previous section as now the number of lines 
entering a vertex v can assume any value s v € N, and if v ^ Vq{9) the corresponding node factor is 

F v = -— } g Sv , (7.4) 

which is bounded proportionally to some constant G to the power s v . Since J2vev(6)( s v ~ 1) = 
\E{9)\ — 1 < k — 1 (by Lemma 3) this produces an overall constant G 2k in the tree value. Also the 
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study of the convergence of both the formal series and the resummed series can then be performed as 
in the previous case, and no further difficulty arises. The constant c appearing after (6.14) becomes 
— g'(co), instead of — 2c , so that still one has c ^ by the assumption (1.6). 



Acknowledgments. We are indebted to Ugo Bessi for many enlightening discussions. 
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